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ABSTRACT 

Humidity  spectra  have  been  measured  with  the  Lyman - 
alpha  humidiometer,  together  with  mean  profiles  of  virtual 
potential  temperature,  wind,  and  humidity  in  the  open 

ocean  environment.   Empirically  derived  expressions  des- 

2 
cribing  the  temperature  structure  parameter,  C„  ,  were 

extended  by  similarity  arguments  to  the  humidity-structure 

2        2 
parameter,  C   ,  and  C   was  related  to  the  stability  para- 
q  q 

meter  Ri.   Using  the  above  measured  parameters,  vertical 
humidity  flux  was  computed  in  two  different  manners ,  and 

a  comparison  was  made.   In  general,  there  was  little 

2 
correlation  between  the  spectrally  analyzed  C   values 

2 
and  Ri.   Results  for  C   can  essentially  be  regarded  as 

q 

a  function  of  z  and  humidity  gradient.   Non-dimensional 

2 
C   results  were  generally  an  order  of  magnitude  smaller 

than  expected.   No  correlation  between  the  two  methods  of 

calculation  of  humidity  flux  was  found. 
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I.   INTRODUCTION 

In  recent  years  the  turbulent  structure  of  the  atmos- 
pheric boundary  layer  has  become  very  important  to  both 
electro -magnetic  and  optical  wave  propagation  studies. 
Small  scale  temperature  and  humidity  fluctuations  are 
important  in  determining  such  phenomena  as  the  intensity, 
fluctuation,  beam  wander,  beam  spread,  and  absorption  of 
a  transmitted  wave.   It  is  also  recognized  that  these 
fluctuations  are  important  in  the  understanding  of  electro- 
magnetic waves . 

Initial  experimental  efforts  to  verify  turbulence  theory 
predictions  were  conducted  over  land.   Measurements  of  wind 
speed,  temperature  and  humidity  fluctuations  in  those  in- 
vestigations were  taken  from  a  stable  platform  with  instru- 
mentation not  encountering  as  severe  conditions  as  were  these 
The  Navy,  however,  is  interested  in  the  prediction  of  these 
parameters  in  the  marine  environment.   Experiments  at  the 
Naval  Postgraduate  School  have  been  providing  the  basis  for 
predictions  of  the  properties  of  the  atmosphere  for  rele- 
vant wavelengths.   The  turbulent  fluctuation  of  temperature 
is  that  parameter  which  is  important  in  optical  wave  propa- 
gation.  Numerous  studies  of  the  turbulent  temperature 
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fluctuations  have  yielded  significant  results.   In  this 
analysis,  another  facet  of  the  atmospheric  boundary  layer 
turbulent  regime  is  considered,  namely  humidity. 

It  was  the  purpose  of  this  thesis  to  examine  the  turbu- 
lence properties  of  water  vapor  in  the  atmosphere.   It  was 
examined  in  a  fashion  similar  to  that  used  for  temperature, 
since  both  parameters  are  passive  scalar  contaminants  in 

the  atmosphere.   An  analysis  of  the  humidity- structure 

2  2 

function,  C   ,  was  made,  and  an  attempt  to  relate  C   to  a 
q  q 

stability  parameter,  the  Richardson  number  was  made.   Like- 
wise, the  turbulent  transfer  of  the  vertical  humidity  flux 


(w'q')  was  examined. 

This  study  was  based  on  data  obtained  during  April  1976 
under  open  ocean  conditions  aboard  the  R/V  Acania.   The 
data  were  examined  on  the  basis  of  expressions  used 
previously  by  Wyngaard  et  al.  (1971)  and  Friehe  et  al, 
(1975). 


13 


II.   THEORETICAL  BACKGROUND 

A0   GENERAL 

Based  on  analyses  of  Corrsin  (1951)  the  variance  spectra 
of  a  passive  scalar  property,  such  as  humidity,  at  large 
Reynolds  numbers  should  have  an  inert ial  subrange  of  the 

form 

-1/3   -5/3 
S(k)  =  /3X     e    i/J  k  D/J   ,  (1) 

q        q 

where  y5  is  an  empirical  constant  evaluated  to  be  0.25. 

The  temperature  spectrum  in  terms  of  the  structure  function 

2 
parameter  C   ,  has  the  form 

ST(k)  =  0.25  CT2  k"5/3   .  (2) 

Similarity  arguments  lead  to  a  parallel  expression  for  the 

2 
humidity  structure  function  parameter,  C 

2   -5/3 
S  (k)  =  0.25  C   k  '    .  (3) 

q        q 

Equation  (3)  defines  the  one -dimensional  humidity  spectrum 
which  by  definition  is  the  Fourier  transform  of  the  corre- 
lation function  with  separation  r  in  the  streamwise 
direction. 
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B.   CALCULATION  OF  THE  HUMIDITY -STRUCTURE  PARAMETER 

2 
Equation  (3)  can  be  solved  for  C   as 

q 

2  5/3 

Cq  =  4  S  (k)  kD/J  (4) 

where  k  is  the  streamwise  component  of  wave  number.   Since 
the  humidity  fluctuations  are  measured  at  a  fixed  point  in 

the  flow,  the  resultant  spectra  are  realized  at  a  temporal 

2 
frequency,  f.   To  obtain  C   ,  the  temporal  (f)  and  space 

(k)  scales  are  assumed  to  be  related  by  Taylor's  frozen- 
field  hypothesis, 

k  =  2  7rf/U  (5) 

where  U  is  the  mean  wind  at  the  measurement  level.   "Frozen" 
implies  that  the  turbulence  pattern  remains  unchanged  as  it 
sweeps  past  the  probe.   This  yields 

C  2  =  4  S  (kH^-)5/3  f5/3.  (6) 

3     q    u 

This  is  the  expression  that  was  used  to  compute  the  humidity- 
structure  parameter  from  humidity  variance  spectra,  S  (k) . 

C0   STABILITY  CONSIDERATIONS 

Present  boundary  layer  turbulence  theory  has  its  basis 
on  empirical  results  obtained  by  Monin  and  Obukhov  (1954). 
The  Monin-Obukhov  theory  defines  a  scaling  length,  L, 
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proportional  to  the  level  where  mechanical  and  thermal 
production  of  turbulent  kinetic  energy  are  equal.  The 
Monin-Obukhov  length  is  defined  as 


L  =   q*T„  <7) 


where  T^.  -  T*  +  0.61  q^  T,  q^  =  -w'q'/U*,  and  K   =  0.35 
(von  Karman's  constant).   The  ratio  of  the  height  of  the 
measurement  to  the  Monin-Obukhov  length,  z/L,  serves  as  a 
stability  index. 

Recent  observational  experiments  by  Businger  et  al 
(1971),  as  well  as  others,  have  provided  a  relationship 
between  the  Richardson  number,  Ri, 

Ri  =  »f  nV/^  (8) 

9(  du/  dz)  ^ 

and  the  Monin-Obukhov  length,  L,  where  0v  is  the  virtual 
potential  temperature.   Figure  (1)  from  Businger  et  al 
(1971)  illustrates  this  determination  quite  well.   The 
latter  results  substantiated  expressions  for  the  relation- 
ships between  z/L  and  Ri  proposed  by  Dyer  and  Dicks  (1970) 
for  unstable  conditions.   These  expressions  are,  respectively, 
z/L  =  Ri  (9) 
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Here  a    is  an  empirically  derived  constant  equal  to  0.5. 

It  is  important  to  note  that  in  these  expressions  z/L 
approaches  <*  as  the  Richardson  Number  approaches  a  critical 
value  of  0.21.   This  suggests  that  as  stability  increases, 
the  flow  becomes  essentially  non-turbulent  and  the  effect 
of  mechanical  turbulence  becomes  negligible. 
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Figure  1. 
The  dependence  of  Richardson  Number  on  stability. 
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Applying  the  definition  of  Eqn  (4)  to  the  general 
expression  of  Eqn  (1),  it  is  seen  that 

Cq2=  ^  V"1/3  (U) 

2 

C   as  defined  by  Eqn  (11)  enables  indirect  estimates  of 

2 

C   to  be  obtained  based  on  mean  conditions.   This  is  because 

q 

€  and  X     are  easily  related  to  boundary  fluxes  and  profiles, 

if  steady  and  horizontally  homogeneous  conditions  exist. 

2 
Expressions  which  relate  C   to  mean  properties  of  the  bound- 
ary layer,  such  as  Ri  or  z/L,  are  desirable  because  the  small 
scale  measurements  are  impractical  to  obtain  in  most  opera- 
tional or  tactical  situations. 

2 
Expressions  relating  C   to  €  and  X    ,  and  similarity 

theory  predictions  for  the  dependence  of  €     and  X     on 

In 

momentum  and  humidity  fluxes  were  obtained  from  similarity 
arguments  by  Wyngaard  et  al  (1971) : 


TV     ■   fi(z'L>  <12> 


*q 


q*u*z     2 


=  f0(z/L)  (13) 


k 


where  U*=(-u'w')2,  q*=-w'q7U*,  L=-TQU^  /kgw'T'  and  f^z/L) 
and  fo(z/L)  are  empirically  determined  functions. 
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Direct  substitution  of  Equations  (12)  and  (13)  into 
Equation  (1)  yields  an  expression  of  the  form 

2     2   -2/3 


C 


q*   z  *'J  f3(z/L)  .  (14) 


f^(z/L)    is   defined  by   the  combination  of   f   (z/L)    and  f?(z/L) , 
and   is   given  by 

(4.9    [1    -   7(z/L)]    "2/3;    z/L<0 
f3(z/L)    =  (15) 

(4.9    [1+2. 8 (z/L)]  ;    z/L  >  0 

Furthermore,  since  z/L  and  Ri  are  functionally  related, 
as  described  previously,  a  parallel  dependence  of  Ri  can 
be  obtained,  which  is  (Wyngaard  et  al) 

Cq2    =   z4/3(  aq/az)2    f4(Ri)  (16) 

This  final  expression  provides  a  desired  dependence  of 

2  _ 

C   on  more  readily  measured  parameters  (z,  dq/  3 z) ,  and 

Ri). 

Results  for  the  non-dimensional  temperature  structure 
function  parameter  derived  over  open  ocean  conditions  by 
Hughes  (1976)  and  its  comparison  with  the  Wyngaard  et  al 
prediction  curve  are  shown  in  Figure  2. 
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CALCULATION  OF  VERTICAL  HUMIDITY  FLUX 

It  is  possible  to  determine  the  vertical  humidity  flux 


(w'q')  from  measurements  of  wind,  C   ,  and  vertical  humidity 
gradient . 

1.   The  Spectral  Method 

The  spectral  method  yields  the  following  relation- 
ship for  the  viscous  diffusion  of  humidity,  X    , 


X      =  (w'q')  bq/   dz  (17) 

q 

Here,  the  humidity  variance  balance  is  determined  for  steady 
state,  horizontally  homogeneous  conditions,  and  q  is  the 
specific  humidity. 

Combining  Equations  (11)  and  (17)  and  solving  for 


(w'q')    yields 


2        1/3 

^  ■  hrrwr-z  (18) 


Eqn  (18)  needs  to  be  simplified  to  use  with  the  measured 
data.  This  simplification  will  be  described  in  the  next 
section. 

2.   Direct  Profile  Method 

The  acknowledged  expression  for  the  gradient  of 
humidity  pertinent  to  a  steady-state,  horizontally  homo- 
geneous atmosphere  can  be  written  as: 


21 


dq/  dz  =  q*/  Kz  (19) 

Applying  the  definition  of  the  scaling  humidity,  namely 


q*  -  ^iLSLL  (20) 


to  Eqn  (19)  results  in 


>*'•»« -fit  <2« 


When  expressed  in  finite  difference  form,  a  simple  approxi- 
mation to  ^q/^z  is 

AS.I      * £3.   -  q(22)-q(zl)  (22) 

a  z  lz~   4  z   z«1   (z^/z..) 

where  the  computed  number  applies  to  the  height  z^.   Apply- 
ing this  finite  difference  approximation  to  Eqn  (21)  and 


solving  for  (-w'q1)  yields 


-w^v  =  *u*(q(z2>-q(zi>>  (23) 

In  (z2/z1) 


This  is  the  expression  that  was  used  to  compute  vertical 
humidity  flux  from  the  vertical  profile  of  mean  specific 
humidity. 
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III.   THE  EXPERIMENT 

A.   THE  PLATFORM  AND  LOCATION 

Observations  were  made  aboard  the  R/V  ACANIA  anchored 
off  Monterey,  California  in  Monterey  Bay.  Measurements  were 
made  at  multiple  levels  on  two  masts  spatially  separated  on 
the  forward  deck  of  the  ship.   An  illustration  of  the  sensor 
locations  of  these  masts  is  shown  in  Figure  (3) .   The  vane 
and  probe  arrangement  appears  in  Figure  (4) .   The  vane  main- 
tained a  one-dimensional  profile  of  the  wind  during 
measurements . 

The  location  of  Monterey  Bay  provided  an  ideal  site  for 
the  experimental  program  as  shown  in  Figure  (5) .   Open  ocean 
differs  from  land  in  the  effects  of  wave  action  on  turbu- 
lence, in  the  nature  of  the  aerosols  and  fog,  and  in  the 
humidity  fluctuations.   These  conditions  can,  of  course,  be 
best  obtained  far  at  sea.   The  cost  of  such  activities  makes 
it  desirable  to  work  near  land.   Pt.  Pinos  and  Monterey  Bay 
provide  very  nearly  the  ideal  situation. 

Pt.  Pinos  projects  northwest  from  the  mainland  toward 
the  prevailing  northwest  wind.   Even  under  storm  wind  condi- 
tions, the  wind  comes  from  the  southwest,  still  bringing 
sea  air  toward  shore.   The  ocean  depth  surrounding  Pt.  Pinos 
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Figure  5. 
Experimental  site  in  Monterey  Bay 
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is  shallow  enough  that  the  ACANIA  can  anchor  at  any  range 
up  to  about  15  km  to  the  north  and  northeast. 

B.   INSTRUMENTATION 

1.   Wind-Temperature  Measurement  Systems 

The  mean  wind  measurements  were  made  with  a 
Thornthwaite  Associates  cup  anemometer  wind  profile  register 
system,  model  number  104.   In  operation  the  shaft  of  a  three 
cup  anemometer  unit  serves  as  the  shutter  between  a  light 
source  and  photocell  for  each  revolution.   The  cups  are 
plastic  cones  reinforced  with  aluminum  frames.   They  are 
attached  to  the  rotating  shaft  by  stainless  steel  tubing 
spaced  at  120  degree  intervals  about  the  shaft  as  shown  in 
Figure  (6) .   The  three  cup  assembly  sets  along  with  the 
other  sensors  were  positioned  at  four  levels  on  the  bow 
with  electrical  leads  to  the  after  deck  house  laboratory. 
The  sets  have  the  characteristics  of  low  starting  speeds 
with  a  small  amount  of  internal  friction  which  aids  in 
checking  inert ial  overshoot. 

Temperature  sensitive  quartz  crystal  probes,  (Hewlett 
Packard  model  HP -2850)  were  used  to  measure  mean  temperatures 
at  each  level.   RF  signals  from  the  crystal  probes  and  from 
a  reference  oscillator  were  mixed  in  the  HP-2801A  readout 
unit  to  produce  a  beat  frequency  whose  signature  can  be 
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Figure  6. 
C.  W.  Thornthwaite  Anemometer  Cups 
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analyzed  to  within  0.001  degrees  centigrade  per  hertz.   Each 
sensor  simultaneously  received  pre-experiment  calibration 
against  a  platinum  resistance  wire  thermometer  in  a  tempera- 
ture controlled  circulating  water  bath  over  the  expected 
temperature  range.   The  accuracy  in  achieving  a  0.005  degree 
centigrade  correction  factor  was  a  constant  for  each  probe. 
A  3.7  meter  flexible  coaxial  cable  is  permanently  attached 
to  the  sensor  head  and  the  mast  mounted  probes  are  housed 
in  an  aspirated  shelter  as  depicted  in  Figure  (7) .   Tempera- 
ture values  were  automatically  recorded  on  a  printer  tape. 

The  temperature  fluctuations  were  measured  using  a 
bridge  developed  by  personnel  at  GTE  Sylvania,  the  GTE 
Sylvania  Model  140.   The  system  was  slightly  modified  for 
use  in  this  study,  and  is  described  by  Karch  (1976). 
2.   Humidity  Measurement  Systems 

To  determine  their  accuracy  and  reproducibility, 
detailed  examination  in  laboratory  conditions  has  been  made 
of  the  calibration  characteristics  of  two  humidity  sensors, 
the  Hygrodynamics  Digital  Hygrosensor  and  the  Electro- 
Magnetic  Research  Corporation  Lyman-alpha  Sensor.   These 
calibrations  were  based  on  ambient  humidity  variations 
(from  24%  RH  to  507o  RH  over  a  one  month  period)  ,  as  measured 
with  a  Bendix  Psychron  Wet-Dry  Bulb  Psychrometer,  and  con- 
trolled humidity  (from  40%  RH  to  99%  RH)  as  measured  by  the 
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C.C.  Breidert  Company 
Air-X-Houster  Type  6L 


Dunmore-type  Lithium 
Chloride  Sensor 


Quartz  Thermometer  Probe 

Figure  7. 
Photographs  of  aspirated  shelter,  hygrometer 
and  thermometer  probe 
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wet-dry  bulb  thermocouple  in  a  2  x  3  x  3  ft.  fog  chamber 
adapted  for  these  experiments.   The  air  circulated  through 
the  chamber  can  be  pre-humidified  in  a  column  of  wet  glass 
beads  giving  nominal  humidity  of  about  807=,.   Higher  humidi- 
ties (up  to  99%  RH)  can  be  obtained  by  introducing  varying 
amounts  of  water  in  the  form  of  a  fine  mist  into  the  chamber. 
In  general,  the  humidity  in  the  chamber  can  be  stabilized  for 
periods  of  several  hours  and  is  uniform  throughout  the 
chamber  to  within  1%   RH. 

a.   Hygrodynamics  Digital  Hygrosensor 

The  hygrosensor  is  a  Dunmore-type  lithium 
chloride  sensor  whose  resistance  varies  in  proportion  to 
the  relative  humidity  to  which  it  is  exposed.   It  is  a  slow- 
response  sensor  and  is  generally  utilized  in  mean  system 
measurements.   The  standard  procedure  of  calibrating  these 
sensors  in  a  small  chamber  with  a  saturated  salt  solution 
of  known  vapor  pressure  was  not  satisfactory  for  a  program 
of  calibration  of  eight  sensors.   This  dissatisfaction  was 
primarily  due  to  a  lack  of  confidence  in  the  ability  to 
reproduce  humidity  conditions  in  consecutive  sensor  calibra- 
tions.  Attempts  to  calibrate  all  eight  sensors  simultaneously 
in  a  large  volume  chamber  were  complicated  by  temperature 
drifts  that  made  it  difficult  to  assure  equilibrium  satura- 
tions necessary  to  the  accuracy  of  the  technique. 
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As  a  result  of  the  above  problems  it  was  decided 
to  use  the  humidity  controlled  chamber  previously  described. 
The  quantitative  results  of  their  calibration  are  given  in 
Table  I;  in  general,  the  eight  sensor  average  agrees  with 
the  psychrometer  standard  to  -0.4  +  2.97o  RH  for  measurements 
taken  with  equilibrium  times  of  order  one  hour.   The  indi- 
vidual total  humidity  accuracy  of  the  sensor  is  of  less 
importance  that  the  relative  consistency  within  the  group 
of  sensors,  since  we  use  them  in  a  four-level  system  to 
determine  the  humidity  gradients.   In  this  respect,  they  can 
be  calibrated  to  +  1%  RH,  subject  to  assurances  of  sufficient 
time  for  response. 

The  suitability  of  the  hygrosensors  for  ship- 
board measurement  of  humidity  gradients  is  reduced  by  two 
inherent  properties  of  the  device.  An  aspirated,  teflon- 
coated  sensor  has  a  time  constant  of  about  30  minutes,  but 
this  varies  from  sensor  to  sensor.   Therefore,  the  humidity 
difference  read  between  two  sensors  will  be  unreliable  if 
the  humidity  is  changing  faster  than  about  10%  RH  per  half- 
hour.  A  more  serious  deficiency  appears  if  the  sensors  are 
exposed  to  humidities  about  957o  RH.   At  those  high  humidi- 
ties, the  sensors  become  sluggish  and  exposure  to  99%  RH 
may  cause  the  sensors  to  be  useless  for  as  long  as  several 
days. 


32 


TABLE  I 
Hygro dynamics  Digital  II  Calibration 


Wet- 
Dry 
Bulb 
%RH 

Ave 

Hygro 

%RH 

Deviation  of 
12            3 

Sensor 
4 

from 
5 

Ave  Hygro  7> 
6             7 

RH 
•   8 

33. 

3.16 

1.4 

-1.4 

-0.6 

1.5 

-0.2 

-1.6 

-0.9 

1.8 

36. 

37.2 

-0.5 

-0.8 

-0.3 

1.6 

0.5 

-1.4 

-0.1 

1.2 

44. 

41.2 

-0.8 

-0.4 

-0.3 

0.6 

0.9 

0.4 

-0.1 

0.3 

NA 

41.5 

-0.5 

-0.4 

-0.3 

0.5 

0.9 

-0.4 

-0.1 

-0.1 

NA 

42.6 

-1.2 

-0.5 

-0.5 

-0.9 

1.8 

-0.4 

-0.2 

0.1 

NA 

65.8 

-0.2 

-1.4 

-2.5 

0.8 

1.8 

-0.3 

1.3 

0.3 

74.8 

73.6 

0.2 

0.2 

-1.3 

0.5 

-0.2 

0.4 

-0.4 

0.8 

88.3 

85.0 

-0.7 

0.9 

-0.4 

0.1 

-0.6 

0.0 

0.4 

0.2 

82.6 

86.7 

-1.4 

-0.6 

-0.5 

0.1 

-0.2 

0.3 

0.7 

0.3 

82.6 

87.0 

-1.0 

0.5 

-0.7 

0.0 

0.3 

0.7 

0.0 

0.2 

85.0 

88.0 

-0.4 

1.1 

-0.1 

0.1 

-1.8 

1.1 

-0.6 

-0.1 

93.9 

89.6 

0.1 

0.8 

-0.2 

0.3 

-2.2 

2.1 

-1.5 

0.4 

94.4 

92.3 

-0.6 

1.1 

-0.5 

-0.3 

-0.4 

1.2 

-0.5 

0.4 

94.4 

93.6 

-0.8 

1.0 

-0.8 

-0.3 

0.4 

-0.2 

0.8 

-0.4 

Ave  Deviat  -0.5   -0.0   -0.6   0.3   0.1   0.1   -0.1   0.4 
0.7   0.9   0.6   0.7   1.1   1.0   0.7   0.6 


The  average  deviation  of  the  wet -dry  bulb  from  the  Ave  Hygro 
reading  is  0.47,  RH  with  a  standard  deviation  of  2.97,  RH. 
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b.   Electro -Magnetic  Research  Corporation 
Lyman -alpha  Sensor 

The  Lyman-alpha  sensor  measures  humidity  as  a 

function  of  the  absorption  of  ultraviolet  light  by  the  Model 

o 
1215A  Lyman-alpha  transition  of  hydrogen  in  water  vapor.   It 

consists  simply  of  a  UV  source  tube  and  detector,  separated 

by  an  absorbing  air  gap. 

The  detector  voltage  is  proportional  to  the 

s 

total  light  transmitted  and  is  related  to  the  absolute 
humidity,  q  (in  mb) ,  by  the  equation 

V  =  VQ  exp  (-Aq)  •        (24) 

The  Lyman-alpha  is  a  fast  response  device  and  is  specifically 
used  for  measurement  of  humidity  fluctuations,  q'.   Eqn  (24) 
can  be  manipulated  to  relate  q'  to  voltage  fluctuations  V, 
giving 

q'  -  --LJC  (25) 


The  quantity  V   in  Eqn  (24)  is  the  sensor  output  voltage  in 
o 

dry  air  which  can  be  found  by  placing  the  detector  in  de- 
humidified air. 

Because  water  soluble  windows  (LiF  or  MgF  )  are 
required  for  transmission  in  the  ultraviolet,  V   is  a  highly 
variable  quantity,  a  property  which  makes  the  Lyman-alpha  a 
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poor  device  for  measurement  of  total  humidity,  q.   However, 
Eqn  (25)  shows  that  the  fluctuating  part,  q' ,  depends  only 
on  k  ,  which  is  a  property  of  water  vapor  and  the  gap  setting 
only  and  is  independent  of  the  window  condition  so  that  q' 
can  be  reliably  measured.   The  Lyman-alpha  has  been  cali- 
brated in  dry  air,  ambient  air,  and  in  a  humidity  chamber 
and  it  was  found  that  A.  =  .225  mb"1  for  a  source  current  = 
50  fia   and  a  1  cm  air  gap.   The  calibration  curve  appears 
in  Figure  (8) . 

The  feasibility  of  the  Lyman-alpha  sensor  for 
shipboard  humidity  fluctuation  measurements  is  subject  to 
considerations  of  the  survivability  of  the  windows  in  the 
ocean  environment  of  high  humidity,  sea  spray,  and  rain. 
In  the  laboratory  humidity  chamber,  it  was  observed  that 
VQ  would  decrease  by  a  factor  of  four  after  a  three-day 
exposure  to  85%  humidity.   Consequently,  it  will  be  neces- 
sary to  provide  better  physical  protection  for  the  sensors 
during  routine  operations. 
3.   Data  Recording 

All  mean  data  recorded  for  this  study  was  logged 
using  an  NPS  developed  micro-processor  based  data  acquisi- 
tion system,  the  MIDAS  (Microprogrammable  Integrated  Data 
Acquisition  System)  .   It  is  fully  described  by  Atkinson 
(1976) .   In  this  study  a  ten  minute  averaging  period  was 
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Figure  8. 
Calibration  graph  of  Lyman-alpha  humidiometer 
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used  to  define  the  mean  data.   Temperature,  relative  humidity, 
and  mean  wind  data  was  punched  on  paper  tape  by  the  MIDAS 
unit;  the  data  was  later  transferred  to  computer  cards  for 
further  analysis. 

Voltage  signals  corresponding  to  measured  fluctuat- 
ing humidity,  q',  were  recorded  on  magnetic  tape  using  a 
Sanborn  Model  3950  fourteen  channel  tape  recorder.   In  addi- 
tion, the  average  voltage  output  of  the  Lyman-alpha  sensor 
was  recorded  on  a  brush  type  chart  recorder.   The  latter 
was  necessary  to  define  V  in  the  voltage  to  humidity  con- 
version expression,  Eqn  (25) . 

C.   ANALYSIS  PROCEDURES 

Data  analysis  was  carried  out  for  data  collected  during 
a  four  day  cruise  in  April  1976.   Mean  data  was  obtained  for 
each  of  the  four  levels  during  the  entire  cruise.   However, 
only  two  Lyman-alpha  humidiometers  were  available  so  q1  data 
was  collected  at  only  two  levels.   These  two  instruments 
were  placed  at  levels  6.6  and  13.9  meters  above  mean  water 
level.   Due  to  the  failure  of  a  source  tube  in  the  instru- 
ment at  the  upper  level,  q'  data  was  available  from  only 
the  6.6  meter  level.   Table  II  summarizes  periods  examined 
and  data  available  for  each  period. 
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TABLE  II 
Summary  of  Data  Periods 


Date 

Time 

Number  of 
Levels 
q/q' 

Number  of 
Profiles /Spectra 
q/q' 

1976 

27  April 

1446-1930 

4/1 

12/6 

28  April 

1320-1939 

4/1 

33/24 

29  April 

1410-2240 

4/1 

42/24 

30  April 

1033-1138 

4/1 

21/13 

All  data  were  first  edited  for  gross  errors  or  incon- 
sistencies due  to  equipment  malfunctions.   The  criterion  at 

this  point  for  retaining  or  rejecting  data  periods  depended 

2 
on  whether  Cq  values  could  be  computed  from  humidity  spec- 
tra for  the  time  periods  involved. 

Lyman-alpha  data  were  analyzed  on  the  basis  of  the  strip 
charts  to  give  an  average  voltage  output  at  ten  minute  inter- 
vals.  Relative  humidity  data  from  the  MIDAS  output  was 
converted  to  specific  humidity  q,  by  use  of  the  governing 
meteorlogical  equations  for  moist  air.   Virtual  potential 
temperatures,  9v,  were  calculated  using  IBM  360  routines. 
The  values  of  q   and  0v  were  then  plotted  versus  the  logarithm 
of  height.   Best  fit  lines  were  drawn  to  the  points  from 
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which  slopes  were  picked  off  and  applied  to  the  expressions 
developed  by  Wyngaard  et  al. 

1.  Profile  Analyses 

Virtual  potential  temperature,  9v,  and  specific 
humidity,  q",  values  were  plotted  on  semi-logarithmic  paper, 
since  0v,  and  "q  are  parameters  which  vary  logarithmically 
with  height.   A  best  fit  straight  line  was  then  drawn  to 
the  data  points.   The  procedures,  of  course,  were  subjective, 
and  in  many  instances  several  different  slopes  could  be 
obtained  from  just  one  graph.   Hence  a  criterion  used  was 
not  to  give  single  data  points  too  much  weight  in  determining 
the  best  fit  line.   Consequently,  the  line  drawn  represented 
a  most  probable  position  between  data  points  as  illustrated 
in  Figure  (9) . 

Difficulty  arose  whenever  anomalous  points  appeared 
in  a  graph  as  in  Figure  (10)  where  the  specific  humidity  at 
the  second  level  was  obviously  inconsistent  with  the  other 
three  levels.   In  a  case  such  as  this,  the  inconsistent 

point  was  merely  ignored. 

2 

2.  C_  Analysis  Procedures 

Analog  spectral  analyses  procedures  were  applied  to 
the  humidity  fluctuation  data  to  obtain  variance  spectrum. 
Appropriate  21  minute  segments  of  data  recorded  on  magnetic 
tape  were  recorded  into  an  EMR-Schlumberger  Model  1510 
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Figure  9. 
Typical  profiles  for 

a)  specific  humidity  -  30  April  76  (1210) 

b)  virtual  potential  temperature  -  29  April  76 
(1052) 
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Figure  10. 
Profile  of  27  April  76  (1817)  showing  an 
anomalous  specific  humidity  reading  at  the 
6. 6  meter  level. 
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Digital  Spectrum  Analyzer.   The  procedures  required  to  con- 
vert spectral  values,  obtained  with  the  digital  spectrum 

2 
analyzer,  to  engineering  units  and  then  to  obtain  C   values 

from  appropriate  spectra  are  described  in  the  following 

subsection.  A  parallel  discussion  for  temperature  spectra 

appears  in  Hughes  (1976). 

a.   Scaling  Spectral  Plots 

A  necessary  procedure  was  to  scale  the  spectral 
plots  to  relate  RMS  input  voltages  to  power  spectral  densi- 
ties (PSD) ;  variance  per  unit  frequency.   To  obtain  power 
spectral  density  levels,  corresponding  to  RMS  voltage  inputs, 
calibrated  scale  charts  had  to  be  constructed. 

For  purposes  of  the  x-y  plot  format  of  the 
analyzer  output,  the  RMS  voltages  were  converted  to  y=log-.Q 
(voltage)  units  and  a  graduated  scale  was  constructed  so 
that  the  logarithm  of  volts  RMS  could  be  interpolated  from 
spectral  plots.   The  value  of  the  vertical  scale  (y)  was 
adjusted  for  each  spectrum  as  a  function  of  both  input  gain 

and  spectral  gain.   These  values  were  then  converted  to  PSD 

2 
levels  for  use  in  calculating  C   values. 

The  relation  between  Volts  RMS  to  PSD  units 

used  was 

S(f)  |PSD  units}  =  ^f1:  ^eV!1.!^S)2  X  (10^)2   (26) 
'         '      1.5  Bandwidth 


42 


where  Bandwidth  =  Freq .  Range 

No.  Channels 

(  =  256_Hz  =  l   H  x 
^    256        Z; 


and  cal.  level  Vrms  =  voltage  at  y=0 


(  =  1  Vrms  for  3.16  V  input 
setting) . 


Amplitude  scaling  calibrations  were  accomplished 
using  externally  generated  "white  noise"  signals  of  1  Volt 
RMS  with  a  frequency  range  from  0.1  Hz  to  1000  Hz  (giving  a 
PSD  of  10"3V2/Hz).   Setting  a  3.16  V  (lOdB)  input  on  the  EMR 
1510  Digital  Spectrum  analyzer  insures  that  1  Volt  corre- 
sponds to  y=0.   An  example  of  such  a  calibration  plot  for 
a  frequency  range  of  1-200  Hz  (the  frequency  range  of  inter- 
est) is  shown  in  Figure  (11). 

b.   Computation  of  Turbulence  Parameters  from 
Scaled  Spectra 

2 
Values  of  the  turbulence  parameter  C   were 

obtained  from  the  humidity  variance  spectra  on  the  basis 
of  the  formula  for  the  inertial  subrange  in  wave  number 
space  S(k),  Eqn  (1),  which  predicts  a  -5/3  slope  for  the 
spectra  when  plotted  in  a  log-log  format.   Figures  (12), 
(13),  and  (14)  are  typical  spectra  considered  in  the  analy- 
ses.  Humidity  spectra  often  exhibited  slopes  slightly 
different  than  the  expected  -5/3,  except  in  the  frequency 
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Figure  11. 
Calibration  plot  for  spectrum  analyzer 
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Figure  12. 
Typical  humidity  spectra  showing  -5/3 
slope  in  the  range  from  1  Hz  to  10  Hz 
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Figure  13. 
Typical  humidity  spectra  showing  -5/3 
slope  in  the  range  from  1  Hz  to  10  Hz 
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Figure  14. 
Typical  humidity  spectra  showing  -5/3 
slope  in  the  range  from  1  Hz  to  10  Hz 
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range  of  1  Hz  to  10  Hz.   This  feature  of  humidity  spectra 
has  been  observed  by  others,  namely  Friehe  (1975).  A  uni- 
versally accepted  explanation  has  not  been  postulated  to 
account  for  this  phenomena.   Figure  (15)  is  an  example  of 
a  spectra  which  does  not  exhibit  a  -5/3  slope  in  any  fre- 
quency range.  Approximately  10?o  of  spectra  studies  fell 
into  this  category. 

Assuming  -5/3  slopes  for  the  variance  spectra, 
the  intercept  of  the  best  fit  -5/3  slope  with  the  1  Hz 

frequency  line  was  the  spectral  density  denoted  (PSD)  value 

2 
used  to  compute  C   .   The  measured  PSD  value  was  converted 

q 

to  engineering  spectral  density  units  by  the  relation 

S  (f)  =  CT2   •  PSD 
qv  '     L-a 

=  (g/kg/Volt)2  •  Volt2/Hz  =  (g/kg)2/Hz 

(27) 

where  C-r     was  the  calibration  factor  for  the  Lyman-alpha 

humid  iometer . 

Since  the  humidity  fluctuations  were  measured 

at  a  fixed  point  in  the  flow,  the  resultant  spectral  values 

are  defined  at  "temporal"  frequencies  denoted  as  S  (f)  in 

q 

2 
Eqn  (28).   In  order  to  obtain  C   from  Eqn  (6)  temporal  (f) 

and  space  (k)  scales  had  to  be  related.   This  was  accomplished 

by  using  Taylor's  "frozen  turbulence"  hypothesis  discussed 
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Figure  15. 
Humidity  spectra  lacking  a  distinct  -5/3  slope 
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in  section  II. B.   The  following  equation  relates  temporal 
to  wavenumber  variance  spectral  values, 

fSq(f)  =  kSq(k)  =  c2CqV2/3  (28) 

where  Sq(f)  is  the  spectral  density  value  with  units  of 
(g/kg)2/Hz. 

Eqn  (28)  leads  to  the  following  relationship 
between  C2  and  S  (f) , 

cq2  =  4r  fVf)  (29) 


where,  c2  =  0.25,  an  empirical  constant,  and  k  =  2  7rf/U. 

_    2 
Hence,  with  measured  values  of  f,  S -(f) ,  and  U,  Cq  was  then 

determined  from  Eqn  (29)  for  the  level  of  interest,  z. 

3.   Procedures  to  Calculate  Vertical  Humidity  Flux 


Vertical  humidity  flux  (w'q')  was  calculated  using 

two  different  methods ;  one  utilized  the  spectrally  derived 

2 
Cq   values,  and  the  other  utilized  the  profile  of  mean 

humidity.   The  results  were  compared  on  a  linear  plot  to 

determine  if  any  relationship  could  be  observed. 

a.   The  Spectral  Method 


This  method  of  calculating  (w'q')  was  based  on 

2 
C_  values  derived  from  variance  spectra  of  the  Lyman-alpha 

humidiometer  using  Eqn  (18).   In  this  calculation  the  values 
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of  €  which  were  used  were  derived  from  wind  fluctuation 
data  by  Karch  (1976) . 

b.   The  Profile  Method 

This  method  utilized  the  vertical  gradient  of 
mean  humidity  and  the  friction  velocity,  IL,  to  calculate 


(w'q1)  from  Eqn  (23).   The  profiles  of  q"  were  evaluated 
and  values  of  q   were  read  off  at  two  different  levels  from 
the  best  fit  line  drawn  to  all  but  erroneous  points. 

The  values  of  U*  used  in  these  computations 
were  determined  from  the  e  values  obtained  by  Karch  (1976) 
The  final  relation  used  was 

U*  =  (e  K   z)1/3  ,  (30) 

which  appears  to  give  the  most  consistent  estimates  for  U* 
in  the  marine  boundary  layer. 
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IV.   RESULTS 

Observed  humidity  index- structure  parameter  results, 
made  non-dimensional  according  to  the  Wyngaard  et  al  predic- 
tion, Eqn  (16),  appears  in  Figure  (16).   In  the  figure 
individual  data  points  appear  as  dots,  and  averages  over 
Ri  intervals  of  0.25  appear  as  dots  within  a  larger  circle. 
The  error  bars  are  standard  deviations  from  the  mean  within 
each  interval,  while  the  number  at  the  top  of  the  error 
bars  is  the  number  of  observations  defining  the  mean  value. 
For  both  the  stable  (+Ri)  and  unstable  (-Ri)  stratifica- 
tion cases,  there  appears  to  be  little  agreement  and  no 
definite  trends. 

Scatter  in  the  observed  results  in  Figure  (16)  can  be 

o 
attributed  to  scatter  in  both  the  measured  C   values  as 

well  as  3q/  ^z  values.   Deviation  of  humidity  spectra  from 

2 
a  -5/3  slope  also  caused  uncertainty  in  C   estimates. 

The  values  for  d q / d  z  are  important  because  the  humidity 

gradients  were,  in  most  cases,  small,  and  enter  into  the 

2 
normalization  of  C   as  a  squared  value. 

2 
Significantly,  non-dimensional  C   results  obtained  in 

this  study  were  generally  an  order  of  magnitude  smaller 

2 
than  Gj,   results  obtained  by  Wyngaard  et  al  (1971)  and 
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Figure  16 
Results  of  the  spectrally  derived  humidity- structure 
parameter  vs.  the  Richardson  Number 
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and  Hughes  (1976).  All  sensors,  calibrations,  and  govern- 
ing equations  were  rechecked.   All  calculations  were  redone. 
However,  the  final  results  remained  smaller  in  value  than 
expected.   It  was  expected  that  results  would  fall  in  the 
same  order  of  magnitude  as  Cr^   since  both  temperature  and 
humidity  are  scalars.   One  possible  explanation  is  that  the 
value  of  /?  from  Eqn  (1)  is  not  the  same  for  both  tempera- 
ture and  humidity,  as  presently  postulated. 

2 
If  the  scatter  of  the  non-dimensional  Cq  values  is 

ignored,  the  results  can  be  seen  to  be  essentially  constant. 

If  one  is  willing  to  accept  the  error  of  an  order  of  magni- 

_  2 

tude,  then  given  a  gradient  of  q  and  a  height  z,  C   can  be 

computed  according  to  the  expression 

2    4/3  2 

CqZ  =  z  '   (  dq/  dz)Z  (31) 

The  comparison  of  the  vertical  humidity  flux  calculations 


appear  in  Figure  (17).   (w'q')  obtained  from  Cq^  data  are 
denoted  (w'cl,)xJ-a  anc*  vertical  humidity  flux  obtained  from 


profiles  of  dq/  dz  are  denoted  (w'q')p.   There  is  no  ob- 
served correlation  between  the  two  values.   Perhaps  new 
dimensionalizing  parameters  need  to  be  developed  in  order  to 


be  able  to  predict  values  of  (w'q')  accurately. 
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Figure  17. 
Results  of  the  comparison  between  vertical  humidity 
flux  calculated  by  the  profile  method,  (w'q')p,  vs. 
the  vertical  humidity  flux  calculated  using 
C  2  values,  (w'q')L.a 
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Table  III  summarizes  data  used  to  construct  Figure  (16), 
Tables  IV  and  V  contain  values  used  in  the  comparison  shown 
in  Figure  (17).   Values  in  all  three  tables  were  computed 
at  the  6.6  meter  level. 
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TABLE  III 
C  2  and  Ri  Results 


Time 

U 

5q/5z 

Cq2   x  10"3 

Ri 

27  April 

1446-1509 

7.70 

-0.043 

1.48 

-1.26 

1509-1520 

8.17 

-0.045 

0.88 

-1.07 

1520-1540 

7.84 

-0.048 

0.92 

-1.44 

1540-1600 

7.50 

-0.037 

1.08 

-1.27 

1620-1640 

6.25 

-0.045 

0.96 

-2.73 

1910-1930 

4.03 

-0.030 

1.72 

-0.90 

28  April 

1320-1337 

5.16 

2.60 

1337-1355 

5.30 

-0.026 

2.04 

-1.32 

1355-1413 

5.16 

-0.025 

1.72 

-0.93 

1415-1420 

5.06 

-0.028 

1.24 

-1.05 

1433-1451 

5.02 

-0.024 

1.96 

-1.32 

1452-1504 

4.39 

-0.024 

2.08 

-0.95 

1511-1517 

4.88 

-0.023 

2.00 

-1.26 

1526-1533 

4.50 

-0.031 

1.84 

-2.96 

1540-1551 

4.31 

-0.031 

2.60 

-3.20 

1556-1606 

4.40 

-0.024 

2.05 

-3.70 

1613-1621 

5.30 

-0.032 

1.67 

-1.42 

1633-1641 

4.97 

-0.039 

1.52 

-2.19 

1652-1701 

4.73 

-0.038 

1.27 

-0.58 

1720-1738 

4.31 

-0.044 

1.09 

0.19 

1739-1744 

4.03 

-0.051 

1.56 

-0.74 

1757-1807 

4.31 

-0.048 

1.00 

-0.58 

1807-1817 

4.31 

-0.053 

1.29 

-0.58 

1817-1827 

4.08 

-0.057 

1.36 

-0.52 

1827-1837 

4.03 

-0.052 

1.32 

0.24 

1837-1847 

3.91 

-0.052 

1.07 

0.29 

1847-1857 

3.61 

-0.053    • 

1.16 

0.29 

1909-1919 

2.80 

-0.045 

1.72 

1.68 

1919-1929 

2.50 

-0.033 

1.51 

7.01 

1929-1939 

2.50 

-0.028 

1.60 

7.89 

29  April 

1410-1420 

5.88 

-0.020 

2.36 

0.79 

1420-1430 

5.86 

-0.039 

3.36 

0.12 

1442-1452 

5.86 

6.40 

1500-1513 

6.43 

3.45 
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Table  III  Continued 


Time 

U 

5q/5z 

Cq2  x  lO"3 

Ri 

1521-1532 

6.00 

3.10 

1541-1553 

4.50 

3.77 

1600-1610 

5.86 

2.51 

1653-1703 

4.26 

-0.039 

2.36 

-6.15 

1703-1713 

4.92 

-0.024 

1.72 

-4.06 

1713-1723 

4.31 

-0.023 

1.68 

-10.94 

1723-1733 

4.22 

-0.014 

1.71 

-5.17 

1733-1743 

4.26 

-0.011 

2.60 

-4.82 

1830-1848 

2.60 

2.28 

1848-1906 

3.61 

2.41 

1916-1934 

2.71 

2.00 

1934-1952 

2.83 

2.16 

2000-2018 

2.95 

2018-2036 

3.36 

2036-2054 

3.76 

1.27 

2054-2114 

3.81 

1.25 

2125-2145 

3.66 

1.40 

2145-2157 

3.76 

1.60 

2203-2223 

3.76 

1.69 

2221-2240 

4.71 

1.60 

30  April 

1033-1043 

4.73 

-0.047 

2.04 

-1.56 

1043-1053 

4.45 

-0.043 

1.55 

-2.52 

1053-1103 

4.31 

-0.030 

2.44 

-3.91 

1108-1118 

4.40 

-0.043 

1.76 

-2.90 

1118-1128 

4.55 

-0.043 

1.75 

-2.52 

1128-1139 

4.78 

-0.033 

3.00 

0.37 

1148-1206 

5.53 

-0.031 

1.44 

0.74 

1206-1224 

5.53 

-0.036 

1.65 

0.62 

1224-1242 

6.05 

-0.054 

1.04 

0.0 

1242-1250 

4.40 

-0.041 

1.52 

0.56 

1250-1308 

3.26 

-0.034 

1.83 

0.0. 

1310-1328 

3.41 

-0.040 

1.32 

-1.98 

1328-1346 

3.51 

-0.032 

0.83 

-0.96 
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TABLE 
Results, 

0 

IV 

C   2  Method 
10~3)       dq/dz 

w'q' 

Time 

€(X10' 

-3)    V (x 

w'q  ' 

28  Apr  76 

1337-1355 

2.90 

2.04 

-0.026 

0.045 

1355-1413 

3.67 

1.72 

-0.025 

0.042 

1757-1807 

4.55 

1.00 

-0.048 

0.014 

1807-1817 

4.60 

1.28 

-0.053 

0.016 

1817-1827 

4.66 

1.36 

-0.057 

0.016 

1827-1837 

4.76 

1.32 

-0.052 

0.017 

1837-1847 

4.95 

1.16 

-0.053 

0.015 

29  Apr  76 

1653-1703 

4.00 

2.36 

-0.039 

0.038 

1703-1713 

3.50 

1.72 

-0.024 

0.044 

1713-1723 

2.33 

1.68 

-0.023 

0.039 

1723-1733 

2.80 

1.72 

-0.014 

0.069 

1733-1743 

3.04 

2.60 

-0.011 

0.014 

30  Apr  76 

1033-1043 

0.97 

2.04 

-0.047 

0.017 

1224-1242 

2.09 

1.04 

-0.054 

0.010 

1242-1250 

2.09 

1.52 

-0.041 

0.019 

1250-1308 

1.37 

1.84 

-0.034 

0.024 
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TABLE  V 
w'q1  Results,  Profile  Method 


Time  q(z2)-q(z;L) 


u* 

w'q' 

0.19 

0.013 

0.20 

0.014 

0.22 

0.030 

0.22 

0.030 

0.22 

0.030 

0.23 

0.027 

0.21 

0.024 

0.20 

0.009 

0.18 

0.013 

0.19 

0.009 

0.19 

0.009 

0.13 

0.015 

0.17 

0.023 

0.17 

0.020 

0.15 

0.017 

28  Apr  76 

1337-1355  6.15-6.18 

1355-1413  6.27-6.30 

1757-1807  6.77-6.83 

1817-1827  6.78-6.84 

1827-1837  6.79-6.85 

1837-1847  6.82-6.87 

29  Apr  76 

1653-1703  6.30-6.35 

1703-1713  6.41-6.43 

1713-1723  6.39-6.42 

1723-1733  6.33-6.35 

1733-1743  6.40-6.42 

30  Apr  76 

1033-1043  6.62-6.67 

1224-1242  7.11-7.17 

1242-1250  7.12-7.17 

1250-1308  7.08-7.13 
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